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We present ab initio calculations of the bandstructure of graphene and of short zigzag graphene 
nanoribbons by the screened-exchange-LDA method (sX-LDA) within the framework of density 
functional theory (DFT). The inclusion of non-local electron-electron interactions in this approach 
results in a renormalization of the electronic bandstructure and the Fermi velocity compared to calcu- 
lations within local density approximation (LDA) gives good agreement with experiment. Similarly, 
the band gaps in zigzag nanoribbons (ZGNR) are widened by more than 200%, being of similar 
magnitude than bandgaps from past studies based on quasiparticle bandstructures. We found a 
noticeable effect of non-local exchange on the spin-polarization of the electronic ground state of 
ZGNRs, compared to LDA and GGA-PW91 calculations. 



I. INTRODUCTION 



Since the first successful preparation of free-standing 
graphene, an isolated layer of graphite, in 20043, a con- 
siderable amount of both theoretical and experimental 
work have been employed on investigations of its unique 
physical properties. Graphene is a zero-gap semiconduc- 
tor where the tt valence band and the tt* conduction band 
contact at the six Dirac points in the hexagonal graphene 
Brillouin zone. A peculiar property of graphene is the 
linear dispersion E = hvF\k\ {vp is the Fermi velocity) 
of the TT-bands near the K-points, resulting in effectively 
massless Dirac fermions^^. This property, predicted the- 
oretically and confirmed experimentally, is appealing for 
theoretical scientists, as it allows in principle to study 
relativistic Dirac fermions by methods from condensed- 
matter physics. Similarly, the outstanding electroni- 
cal properties turn graphene and its related materials, 
such as carbon nanotubes and graphene nanoribbons, to 
promising materials for application in nanoelectronics. 
For this reason, a thorough understanding of the elec- 
tronic properties of graphene materials are significant. 
A common and efficient method for the study of ground 
state properties in condensed-matter physics are ab ini- 
tio density functional theory calculations within the lo- 
cal density approximation (LDA) or generalized den- 
sity approximation (GGA). Whereas those approaches, 
in all their simplicity, usually yield good results for 
structural properties they routinely underestimate elec- 
tronical properties, e.g. band gaps. Recent re- 
ports of experimental studies of charge carrier dynam- 
ics in graphene by use of angle resolved photoemis- 
sion spectroscopy (ARPES)^, IR measurement^, mag- 
netotransport measurement^, and scanning tunneling 
microscopji^, revealed that common density functional 
theory calculations severely underestimate the Fermi ve- 
locity by 15 -20%. This was attributed to prominent 
many-body effects in graphene such as non-local electron- 
electron and electron-phonon interaction, which are not 
accounted for in the LDA frame and lead to a renor- 
malization of the Fermi velocity. Similarly, calculations 



within LDA underestimate the band gap sizes in both 
armchair and zigzag graphene nanoribbon^ by 50-75%, 
bad performances even for LDA calculations on semicon- 
ductors. Several authors have reported that self-energy 
corrections of the LDA bandstructures using many-body 
GW approximations can amend the shortcomings of 
LDA. For the Fermi velocity, those corrections yield 
values tha t are in good agreement with experimental 
result^nni. However, while quasiparticle corrections pro- 
duce accurate results, the considerable computational ef- 
fort here is an undeniable disadvantage. 
There are several attempts to improve on the short- 
comings of LDA by including non-local expressions 
in the exchange-correlation functional. In case of 
the screened-exchange (sX-LDA) approximation, the 
electron-correlation functional is modeled by a blend of 
LDA and a screened Hartree-Fock exchange potential 
in order to account for electron-electron interacti ons to 
some extend but retaining the advantages of LDAEESI 
It has been shown that the screened-exchange-method is 
successful in describing the band gap s of v arious materi- 
als with accuracy comparable to GW^-'^'^'. 

Motivated by that success, we want to use this paper 
to report calculations on the bandstructure of graphene 
and zigzag graphene nanoribbons employing sX-LDA. 
We found that the calculated renormalization of the 
Fermi velocity in graphene due to electron-electron in- 
teractions compares well with the experim ental value^^ 
and the ones from quasi-particle correction^^ES!^ Further, 
we show that the non-local exchange has a considerable 
effect on the spins-polarization and the electronic band- 
structures of zigzag nanoribbons, which again compare 
well with quasiparticle bandstructure^^. 

II. METHOD 

Our study of the electronic bandstructures of graphene 
and zigzag graphene nanoribbons is based on pseudopo- 
tential density functional theory in the frame of the 
screened-exchange-LDA approximation, LDA and GGA 
in the PW91 form, and was performed by use of the com- 
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TABLE L Ab initio and experimental values of the 
Fermi velocity in the linear valence bands of graphene. 



Method 


Fermi velocity vf (m/s) 


LDA 


0.89 • 10'' 


GGA 


0.88- lO'^ 


sX-LDA 


1.16- lO'' 


GW (CD integration^^ 


1.12 • lO'' 


GW (Random phase approximation^!^ 


1.25 • lO'' 


Experiment (ARPESj^ 


1.0 ±0.05- lO'^ 


Experiment (IRj^l 


1.02 ±0.01 • 10'' 


Experiment (magnetotransport 


1.07 ±0.01 • 10" 


Experiment (STM;F 


1.1 ■ 10" 



putational package CASTEP^. Graphene was modeled 
by the common two-atomic Wigner-Seitz unit cell and pe- 
riodic boundary conditions in the three spacial directions. 
We found that a distance of 6 A between the periodic im- 
ages is sufficient to minimize interlayer interactions and 
represent graphene as a 2D crystal. The action of the 
atomic core were described by standard normconserving 
pseudopotential^, the valence electrons by plane waves 
with a cutoff energy of 750 eV. Integrations in the recipro- 
cal space were done by a Monkhorst-Pack grid^^ of 9x9x1 
equi-distant points in the Brillouin zone of graphene. We 
further optimized the atomic positions until the maxi- 
mum interatomic force were lower than 0.001 eV/A. We 
modeled graphene nanoribbons by a rectangular unit cell 
being continued periodically along the nanoribbon axis 
and used an energy cutoff of 750 eV and a k-point grid of 
1x1x9 points along the Brillouin zone. The cell dimen- 
sions were chosen in a way to minimize interlayer inter- 
action between the periodic images and to maintain the 
calculational effort at the same time. As for graphene, we 
found that a distance of 6 A between periodic images is a 
sufficient compromise. The dangling bonds of the carbon 
atoms at the nanoribbons edges were passivated with hy- 
drogen atoms in order to maintain the sp^-hybridization. 
The geometries of all nanoribbons were fully optimized. 



III. RESULTS AND DISCUSSION 
A. Bandstructure of graphene 

We used these parameters to calculate the electronic 
bandstructures of graphene in the frame of LDA, GGA 
and sX-LDA. 

Figure [l] (a) compares the valence bands of graphene 
from calculations within the local density approximation 
and with screened-exchange. There is a visible renor- 
malization of the bandstructure due to the inclusion of 
non-local correlation effects. Screened-exchange results 
in a general shift of the three valence cr-bands to lower 
energies and that is most significant near the F-point. 



There, the two degenerate bands from sX-LDA have en- 
ergies of -5eV compared to -3.3 eV from LDA. The con- 
duction bands are not affected as much from non-local ex- 
change. This might correlate with studies by Lee et aZ.E^l, 
who argue that the improved band gaps from screened- 
exchange-LDA functionals in the materials they studied 
mainly originate from a marked downshift of the valence 
bands. This lowered energy of the valence bands com- 
pared to the Fermi energy results in a broadened band 
gap by 1.7 eV at the F-point. The 7r-bands from both 
calculated bandstructures coincide fairly well in the K-M 
part of the Brillouin zone. In the F-K and F-M parts of 
the bandstructure, the non-local electron interaction re- 
sults in a visibly increased slope for the sX-LDA 7r-band 
in respect to the 7r-band from LDA. This results in a 
renormalized Fermi velocity 



1 dE 



(1) 



in the linear part of the 7r-bands from sX-LDA compared 
to LDA. The Fermi velocity is a particularly interesting 
variable in semiconductors and conductors, as it is com- 
parable to the velocity of the electrons that contribute 
to electric conduction. Table U shows a number of re- 
cently reported experimental values for the Fermi veloc- 
ity, which are between 1.0 • 10'^ and 1.1 • lO'' m/s. For the 
bandstructure from LDA, we found a Fermi velocity of 
8.9 • 10^ rn/s, which underestimates the reported exper- 
imental Fermi velocities by 11-19% (depending on the 
experiment). Calculations within the generalized gradi- 
ent approximation yield the same result, see Table [l] In 
contrast, the steeper slope in the sX-LDA bandstructure 
results in a Fermi velocity of 1.16- 10^ rn/s, being in good 
agreement with the values from magnetotransport"^ and 
STM'^ measurements. 

Figure [l](b) shows the valence bands from calculations 
within sX-LDA, together with quasi-particle energies at 
selected k-points from GW calculations by Trevisanutto 
et aM The GW calculations build upon LDA wave- 
functions and have been performed by use of a contour- 
deformation (CD) integration, which is known to provide 
the most accurate results in GW quasiparticle correc- 
tions. As with sX-LDA, the correlation effects induce a 
shift of the three cr-bands to lower energies, but the down- 
shift is considerably weaker and of magnitude <1 eV. Un- 
fortunately, the authors did not report any data for the 
conduction band energies. Lee et aZ.^, however, found 
that the improved band gaps from GW result from a 
strong upshift of the conduction bands. The tt band from 
sX-LDA shows a very good agreement with the 7r-band 
in the quasiparticle bandstructure, except for deviations 
near the F-point. The value of the Fermi velocity was 
reported to be 1.12 • 10^ m/s, i.e., only slightly closer 
to experimental values than our result. This encourages 
us to believe that the electronical properties of graphene 
and graphene-based materials can be described in good 
accuracy within the frame of the screened-exchange ap- 
proximation. 
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FIG. 1: (Color online) Comparison of the bandstructure of graphene within the screened-exchange approximation 
with (a) the bandstructure within the local density approximation (LDA) and (b) the quasi-particle bandstructure 
from self-energy corrections within the GW approximation, which was extracted from Ref. [9l 



B. Spin-induced band gaps in zigzag nanoribbons 

The first investigations of the electronical strutures of 
graphene nanoribbons were based on simple zonefolding 
arguments, neglecting effects from the nanoribbon edges. 
In this approximation, one third of all nanoribbons with 
armchair edges and all zigzag nanoribbons should be 
metallic or zero-gap semiconductors^" similar to the 
situation in carbon nanotubes. More elaborate investiga- 
tions employing density functional theory predicted that 
the existence of a large ratio of edge to 'bulk' atoms in 
carbon nanoribbons indeed has a prominent effect on the 
electronical structure. For zigzag nanoribbons, the decid- 
ing factor are edge-localized electron states, which pos- 
sess a large density of states at the Fermi energy. Using a 
Hubbard model, Fujita et a/.^ showed that the resulting 
Fermi instability leads to the formation of an antifer- 
romagnetic ground state in monohydrogenated nanorib- 
bons in contrast to the dielectric graphene. This spin- 
polarized ground state involves the opening of a direct 
bandgap as the system wants to lower the density of 
states near the Fermi level and remove the instabilitj^^H. 

Investigations employing LDA predict that the 
bandgap size is lower than 0.5 eV for all zigzag nanorib- 
bons and scales antiproportionally with the nanoribbon 
width due to the decrease of spin-interaction between the 
atoms at opposite edge^. It might be expected from 
this interaction that nonlocal components of the elec- 
tron exchange are of importance for the bandgap size 
and that LDA, neglecting nonlocal exchange, is not suffi- 
cient to fully describe the electronical properties of zigzag 
nanoribbons. Indeed, our calculations exhibit marked 
changes in the bandgap size due to the inclusion of non- 
local exchange. Figure [2] (a) shows the electron bands of 
the edge states of a small zigzag nanoribbon from LDA 
and sX-LDA calculations. While both bandstructures 
qualitatively exhibit similar dispersion, the bandgap near 
/c=^ is widened by more than 200% from a value of 
0.5 eV in LDA to 1.65 eV in sX-LDA. Our calculated band 



gap energies are comparable to those obtained from GW 
quasiparticle corrections on LDA wavefunctions by Yang 
et al.U^, see Fig. [2] (b), but the quasiparticle bandgaps 
are consistently 0.2 eV larger than our values. Another 
study, employing the B3LYP hybrid potential, reported 
for an 8-ZGNR (16 carbon atoms, w« 16A) a bandgap 
size of 1.34eVf23l, compared to 1.02 eV from sX-LDA and 
1.23 eV from GW. So far, no experimental data exists for 
comparison of those calculated band gaps sizes with real 
values, however, the similar results from different non- 
local approaches suggests that the small values from LDA 
are indeed a severe underestimation of the real bandgap 
energies. There are two different factors for the size of 
the band gap. One is the quantum confinement, which 
induces a inversely proportional width dependence band 
gap energies, the other is the separation of the valence 
and the conduction band by formation of a spin-polarized 
ground state. To estimate the influence of the spins on 
the band gap size, we fitted a function 

to the sX-LDA band gap energies. A is the deviation 
length from the ideal quantum confinement law Eg oc ^ 
and is an indicator for the effect of the spin interaction. 
A = 9.3 A gives good agreement with our calculated val- 
ues. Yang et al. reported a value of A = 16 A , which 
might imply a stronger contribution of the interaction of 
the electrons with the surrounding system in their GW 
calculations. 

We obtained the spin polarization of the electronic 
ground states of our investigated nanoribbons by opti- 
mization of the atomic positions and the electronic spins. 
Starting from a non-spin-polarized geometry, these op- 
timizations resulted in anti-ferromagnetic ground states 
all our calculations within LDA, GGA-PW91 and sX- 
LDA. Figure [s] (a) shows the momenta of the a-spin di- 
rection in case of a 4-ZGNR, a nanoribbon with 8 carbon 
atoms per unit cell. The ground state from LDA calcu- 
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FIG. 2: (Color online) (a) Edge bands of a small zigzag nanoribbon (6 carbon atoms) from calculations employing 
screened-exchange (black lines) and local density approximation (red dots and broken lines) . (b) Comparison of the 
bandgap sizes of zigzag nanoribbons of various widths from sX-LDA (green circles), LDA+GW (red squares) and 
LDA (blue triangles) calculations. The GW results were taken from Ref. [121 



lations exhibits a comparatively weak spin-polarization, 
which is strongly localized at the nanoribbon edges, see 
Fig.[3](b). The calculated spin momentum at the edges is 
0.09 h and rapidly decreases to a almost negligible value 
of 0.01 h at the nanoribbon 'bulk' atoms. This corre- 
sponds to a decrease of 88% over the distance of second- 
nearest neighbours. While the spin-momenta from sX- 
LDA still show a noticeable localization at the edges, the 
spin-polarization over the whole unit cell of the 4-ZGNR 
is more balanced. The spin momentum at the edges 
is 0.21 fi., more than 2.5x the value from LDA, and de- 
creases by 71% into the nanoribbon 'bulk'. As expected, 
GGA improves on the magnetic properties of the mate- 
rial. We found that that the spin momenta from PW91 
are quite in-between the values from LDA and sX-LDA, 
falling from 0.15 h at the edges down to a value of 0.03 h, 
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FIG. 3: (Color online) (a) Comparison of a spin 
momenta of a 4ZGNR (eight atoms per unit cell) from 
calculations employing LDA (red circles), (b) Graphical 
representation of the spin momenta from (a) on the 
nanoribbon atoms. The ratio of the areas of the circles 
representing the spin momenta is the same as the ratio 
between the absolute values of the momenta. 



i.e. a drop of about 80%. The same results hold for 
the momenta of the j3 spins, but from the opposite edge 
and with negative values, and for wider nanoribbons. It 
is evident from the degree of spin-polarization and the 
changed band gap energies that the nonlocal component 
of the electron exchange indeed has a significant effect on 
the electronic properties of zigzag nanoribbons. 



IV. CONCLUSION 

In conclusion, we used the screened-exchange (sX- 
LDA) approximation to calculate the electronic band- 
structure of graphene and zigzag graphene nanoribbons 
and compared the results with the bandstructures from 
LDA and the quasiparticle bandstructures from the GW 
approach. The introduced electron-electron interaction 
in sX-LDA leads to a renormalization of the bandstruc- 
ture and the Fermi velocity in the linear part of the tt- 
band in graphene, as was found in experiments and from 
GW calculations. We report that the Fermi velocity from 
sX-LDA is comparable to the value from GW calculations 
and in good agreement with experimental values. For 
zigzag nanoribbons, our calculations confirmed the signif- 
icance of the nonlocal part of the electron exchange inter- 
action for the degree of spin-polarization and for the size 
of the spin-induced bandgap in the electronic bandstruc- 
ture. The nonlocal exchange results in a larger bandgap 
compared with LDA and comparable to that from GW. 
The spin momenta in the electronic ground state from 
sX-LDA calculations are noticeably higher and more bal- 
anced over the whole unit cell of graphene nanoribbons 
compared to results from LDA and GGA-PW91. As a 
result, we are confident that the screencd-exchange-LDA 
approach is a useful alternative method for the study of 
electronic properties in graphene-relatcd materials. 
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